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1. Introduction. 


In this paper we study the convergence of spectral methods for nonlinear conservation 
laws. Specifically, we consider what is accepted by now as the universal model problem 
for such scalar laws, namely, the inviscid Burgers’ equation 


( 1 . 1 ) 


d d u 2 (x,t) 

- u( x, 1) + 5 I (-L-') = o 


subject to given initial data u (x,£ = 0). Among the basic features of solutions to this 
problem [6], we recall that they may develop spontaneous jump discontinuities (shock 
waves) and hence the class of weak solutions must be admitted; that within this class, 
there are many possible solutions; and that in order to single out the unique ‘physically 
relevant’ solution among them, (1.1) is augmented with an additional entropy condition 
which requires 


( 1 . 2 ) 




The existence of physically relevant shock waves in the solution is reflected by the strict 
(distributional) inequality in (1.2). 

We want to solve the 2?r -periodic problem (1.1), (1.2) by the spectral-Fourier method. 
To this end, we approximate the spectral-Fourier projection of u(x, t), 

N i r* 

(1.3) P*u(x,*) = Y" u{k,t)e ikx , u(k,t) = — / u(x,t)e ikx dx , 

kt^ N 27r J~* 

by an iV-trigonometric polynomial, ujv(x, t), 

N 

(1.4) u N {x,t)= y u k {t)e' kx . 

k=- N 

Starting with 

(1.5) ujv(x,0) = Pnu(x,0) , 

the classical Fourier method [3] lets u^(x, t) evolve at a later time according to the ap- 
proximate model 


( 1 . 6 ) 
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Noting that Pn commutes with differentiation, we can rewrite (1.6) in the equivalent form 

(1.7) iu w (x,i) + ^(i^(x.0) = (/-P«)^(iu5,(x,()) . 


Let us multiply (1.7) by ufj(x,t): since ujy(x,£) is orthogonal to the righthand-side of 

(1.7), we find after integration that 


1 d f w [ w d 1 

(1.8) 2~dt J u 'N{ x it)dx = — J upf(x,t)-^-(—u N (x,t))dx = 


_ U^(X,<) ,X = 7T 

3 '*=-*■ 


Thus /ujy(x,t)dx is conserved in time 

/ IT r IT r IT 

u%(x,t)dx = J u 2 N (x,0)dx < J u 2 (x,0)dx , 
and this yields the existence of a weak limit u(x,t) =tylim Uff(x,t). Does u(x,t) solves 

N—>oo 

our problem? Unfortunately the answer is no. For otherwise, if u(x, t) is a weak solution of 
(1.1), then PpfU 2 N (x,t ) and hence u^(x, ^) should tend weakly to u 2 (x,<), and consequently, 
u(x,t) should be the strong limit of ujv(z,t); but then (1.9) implies that f^ jr u 2 (x,t)dx is 
also conserved in time, and by (1.2) this contradicts the appearance of physically relevant 
shock waves in our solution. 

In practical applications, spectral methods are often augmented with smoothing pro- 
cedures in order to give a helping hand toward their spectral convergence. Indeed, conver- 
gence for smoothed versions of spectral (and in particular pseudospectral) methods, was 
established in the linear case, e.g. [5], [7], [13]. However, arguments similar to the above 
show that with nonlinear problems, convergence of the Fourier method fails despite the 
additional smoothing of its solution. We leave the details for the appendix. Instead, we 
propose here a different way to enforce the convergence of the spectral-Fourier method 
without sacrifying spectral accuracy. This is accomplished by introducing, in Section 2, 
a new type of spectral vanishing viscosity. In Section 3 we prove the convergence of the 
proposed method using compensated compactness arguments, and in Section 4 it is shown 
that the limit solution respect the entropy condition (1.2). In Section 5 we extend our 
discussion to systems of conservation laws, and we show how the spectral vanishing viscos- 
ity can be used to enforce the correct entropy dissipation in such case. Finally, numerical 
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experiments with the proposed method of spectral regularization are presented in Section 

6 . 


2. The Spectral Vanishing Viscosity. 

It is well known [6] that the unique entropy solution of (1.1), (1.2) is the one identified 
with the small viscosity limit of the regularized problem 


( 2 . 1 ) 


^ 7 u «(x,0 + ^(^(x,t)) = e£- 


dt 


dx 


Qfo u '( x ’ *) 


e | 0 . 


With the vanishing viscosity method [9], one replaces the exact derivatives in (2.1) by 
their discrete counterpart, the viscosity coefficient Q is chosen as (a nonlinear) positive 
grid dependent quantity, .and the role of e is played by some fixed, power of the vanishing 
grid size, e ~ (Ax)*, in order to yield an s-order accurate approximation of (1.1). Yet in 
order to respect spectral accuracy, a more delicate viscous regularization is required. To 
this end we consider viscosity coefficients of the form Q = I — P m . The resulting viscosity 
terms are of spectrally small order of magnitude in the sense that for any s > 0 we have 


d_ 

dx 


(I - P m )^u(x,t) 


< t « 

t a (») 11 ll H* + a (») 


m 


in fact these terms are exponentially small in the analytic case [10]. Together with this 
kind of spectral vanishing viscosity, the spectrally accurate Fourier approximation of (1.1) 
amounts to 


( 2 . 2 ) 


± UN ( z ,t) + 


i 1 ~ Pm)-Q^l N {x,t) 


and we raise the question of its convergence as N tends to infinity. Here e = e(V) | 0 
and m = m(N) < N are free parameters which are yet to be determined, subject to the 
spectral accuracy restriction m(N) | oo. In the next two sections we find such admissible 
parameters which provide a positive answer to the convergence question. 
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3. Spectral Convergence to a Weak Solution. 

We consider the approximate Fourier method (2.2), which we rewrite as 


, v d d . 1 2 \ d 

< 31 > at Ufr + = ai 


(r- JV)|»3Sc 


a 


{I ~ Pm) di UN 


= / + //. 


In order to prove convergence of this method we need a couple of a’ priori estimates on its 
solution. To this end, we multiply (3.1) by ujv, 


(3.2) 


312 


h( + £(!«*») = 


3 

a 


= 


a® 




a 

+ £UJV — 

ox 




= in + IV 


and integrate over the 2n period: the integrals of the second and third terms vanish by 
periodicity and orthogonality, and we are left with 


(3.3) 




This gives us the a’ priori bound on the amplitudes of the solution we had before in (1.9), 
and even a little more. More precisely, temporal integration of (3.3) yields 

(3-4) s l r JP - = jlM-.‘ = o)||’, w , 

and hence for u# (x,t) = uk[t)e tkx we have 

(3.5) 


UtfMlI 2 = ^2 |ujfe(t)|^ < Consti , Const\ = ||u (*, t = 0 
t3(l) \k\<N 


2 

t*(x) 


equality (3.4) also gives us the second a’ priori estimate 

(3.6) e||(/- P m )^p-|| 2 =s f ^2 k 2 \u k {t)\ 2 dt < Consti . 

dx ,? oe („.) Jt m< \ k \< N 

Equipped with these estimates we may turn now to the convergence proof of the 
Fourier method (2.2). We will establish spectral convergence for an admissible set of 
parameters e{N) j 0, m(N) f oo, using Tartar’s div-curl lemma [14]. In order to apply the 
latter in our case, we have to verify that the four expressions appearing on the righthand- 
sides of (3.1) and (3.2) are ‘nice’ ones, namely, that these expresssions are the sum of 
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terms, each of which lies either in the compact of H\ Q *(x,t), or - by Murat’s lemma, in a 
bounded set of t). In the following lemmas we collect the necessary estimates in 

this direction. 


We begin with the first term on the right of (3.1). Here, the following estimate whose 
proof is postponed to the end of this section, is essential. 

Lemma 3.1. There exists a constant Constz (depending on Const i and Const2) such 
that form. < £ N we have 


< 3 - 7 > I k 




Next, we use the a’ priori estimate (3.6) to conclude that, as e tends to zero, the 
second term on the right of (3.1) belongs to the compact of H^l(x,t). 

Lemma 3*2 S The following estimate holds 


(3.8) 


dx 


(/ - Pm) 


dux 


dx 


J *,">•«> 


< e\\(I-P m )^r |L (x t) < (e - Const 2 ) 1/2 -> 0 , e | 0 . 


dx Mt <oc ( 

To treat the expressions on the right of (3.2), we first prepare 
Lemma 3 . 3 . There exist constants Constz and Const 23 ( Const 2 3 depending on Const2 
and Constz) such that 

d*u N 


(3.9) 


dx ‘ 


< Constz • m* 


(3.10) - C ° nSt 23 ‘ ( £ V2 + 

Proof. The inequality (3.9) is immediate in view of || 1| ^ < m a '||u 7 v||^ ^ 

and the a’ priori estimate (3.5). To prove (3.10) we invoke the identity 


duff 11 2 
dx *?„<*’♦> 


|| {I ~Pm) 


du N 

dx 


2 


+ p, 


du N 
1 dx 


2 


and use the a’ priori estimate (3.6) to upper bound the first norm, and (3.9) for the second 


one. 
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In the next two lemmas we turn to deal with the righthand-side of (3.2). For its first 
member - which we express as 


(3.11) 


///=«„ — 
OX 


V-p^-A, 


a_ 

dx 


u n{I ~ p n)^n 


- ?£-(i - Pn)\u 2 n = iih + m 2 


we have 

Lemma 3.4. There exists a constant Const 123 (depending on Const 12 and Const 23 ) such 
that the following estimates hold. 


(3.12) 


In ' m k 


u N {I - Pn)^u% 




< Const 123 • |«jv 


!■«<*,«) 


(eN) 


-1/2 


(3.13) \\III2 = -^L{I-Pn)t; u n\\. < Const m • (e 1/2 +m)-(eJV) 1/2 . 

Proof. The first estimate, (3.12), follows from Lemma 3.1 with Const 123 = Const 12 



5 5 


< Const\2 • «at (sN) 


-1/2 . 


together with Lemma 3.3 we also have 

||//Ja|| <||^|| •||(/-P N )ittjS r || 

and the second estimate, (3.13), follows with C123 = Const 1 2 ■ Const 23 - 


Finally, the second member on the right of (3.2), 
IV = «.„(/- p m ) 9 * u " - 

(3.14) 


dx 2 


= e 


dx 


u „(i-p m )^- 


duN, r n \ du N 

- c sr {1 - 


= TV 1 ! + IV 2 , 
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is estimated as follows. 

Lemma 3.5. There exists a constant Const 32 (depending on Const 2 and Const 3 ) such 
that the following estimates hold. 


(3.15) 




tt N (I-Pm) 


dx 


loc 


u,t) 


< Consts2 # un 


, 1/2 




(3.16) ||JV 2 = -e^^-(I- P m )^~-|| < Const 32 {\ + e 1/2 m) . 

OX OX tf ge (s,t) 

Proof. The first inequality, (3.16), follows from the a‘ priori estimate (3.6), 


IV i 


loc 


(*,«) 


<e u N 


>(*.«) 


I-Pm) 


dus 

dx 


L ?ocl"" 


To prove the second inequality we upper bound 


U *IL~(,.«) ’( Const 2 -e ) l/ 2 . 


/V 2 


< e 


loc 


(*»*) 


I-p 

1 * m) dx 


L Ioc (M 




du N 
1 dx 


L loc^> 




and use the a‘ priori estimates (3.6) and (3.9) to obtain (3.16) with Const 32 = Const 2 • 
[Const 3 + 1). 


We are now ready to find the admissible parameters which meet the assumptions of 
the div-curl lemma. By (3.16), the term IV 2 is bounded in L} oe (x,t) if e l / 2 m < Const.; 
choosing e ~ , then by (3.13) the term III 2 is bounded in L lcM provided 

m 2 iV _1 / 2 < Const; choosing m < min (Const. N 1 / 4 , \N) we conclude by Lemma 3.1 
and Lemma 3.2, that the terms I and II belong to the compact of Hfll(x,t). Moreover, if 
us is uniformly bounded, 


(3.17) S Con3t - • 

then by Lemmas 3.4 and 3.5 we have with this choice of parameters that the terms III\ 
and IV\ are also in the compact of H^,l(x,t). This completes our study of the expressions 
on the right of (3.1), (3.2) and the div-curl lemma applies in our case. We summarize by 
stating 
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Theorem 3.0. Consider the spectral approximation (2.2) with parameters (e, m) which 
satisfy 

(3.18) £ = s(JV) ~ Const. N~ 21 * , m = m(N) ~ Const. N 13 , 0 < /? < - . 

4 

Assume that its solution, u s(x,t), remains unifomly bounded, (S.17). Then us(x,t) con- 
verges boundedly a.e. to a weak solution of the conservation law (1.1) 

Proof. Let uO)(x, £) denote the weak limit of u 3 N (x,t ) 

tylim u 3 N (x, t) = «0)(z, t) . 

Applying the div-curl lemma to (3.1), (3.2) we have the relation 

(3.19) = 4^0-^(3T-3(^) 2 , 


which implies strong convergence. To establish this implication, we follow the argument 
of Tartar [15] (see also [14]), who suggests to consider the weak limit of (us — u^ 1 )) 4 , 


w lim (uff — ut 1 )) 4 = «( 4 ) — 4u( 3 ) ut 1 ) + 6u( 2 )(u( 4 )) 2 — 4(u( 1 )) 3 + (ul 1 )) 4 ; 

AT— ►oo 


using (3.19) and rearranging, we find 

tylim (us — u(0) 4 = —3 |u( 2 ) — (ut 1 )) 2 < 0 , 

and hence u( 2 ) = (ut 1 )) 2 . Consequently, us(x,t) converges strongly to u(x,t) = u( 1 )(x,t) 
in Lf oc (x,t), and by (3.1), u(x,t) is a weak solution of (1.1). 


We do not claim that our parametrization (3.18) is optimal. In particular, the re- 
strictive choice of m(N) could be improved as indicated by the numerical tests described 
in Section 6. On the other hand, we note that the e parametrization (3.18) such that 
em 2 ~ Const, yields, in view of Lemma (3.10), 

(3.20a) < Const. , 

11 OX "t? (*,«) 

I oc 

which is in complete agreement with the behavior of the viscous regularization model (2.1), 
where 






< Const. 


(3.206) 
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Our choice of parameters in (3.18) depends heavily on the essential estimate (3.7), 
and we conclude this section with its proof. 

Proof, (of Lemma 3.1.) We should upper bound the norm of 


(3.21) 




JV+P 


= k[ 52 I E «k(0«p-*(0l 2 + \ [ 52 I ^ 

Jt N<p<2N k=p-N Jt -2N<p<-S k=-N 

The first intergral on the right does not exceed 


U £ 


N<p<2N 


p/2 


N 


52+52 |ifi*wi-iv*wi 

K k=p-N k=p/2 


-*/ e 


AT 


N<p<2N 

and using Cauchy-Schwartz inequality this is less than 


E l**(*)HV»MI 

k—p/2 


2 f 52 ( 52 fc2 i«*(*)i 2 - 52 ^i«p-fc(or 

Jt N<p<2N \k>p/2 k>p/2 

According to our assumption m < ^N. Hence for p > N we have p/2 > m, and by the a‘ 
priori estimates (3.5), (3.6), the last expression is bounded from above by 



2- Max ( E 3-EIV‘WI 2 ) ■ / E 

\N<p<2N P k J Jt k>m 


Consti 


<2- 52 4.Afax||«,(,t)|| 2 •Il( / -- P -)^ L IL <T7 -Const! 

P 2 t 11 V " OX "l (*,«) JV 

N<p<2N y ,oc 

The second integral on the right of (3.21) can be treated similarity and Lemma 3.1 follows 

with Const 12 = 4 • (Const i • Const 2 ) 1 / 2 . 
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4. Spectral Convergence to the Entropy Solution. 

In the last sections we have seen that the spectral approximation, u^(z, £), has a 

strong limit, L 2 oe — slim us(x,t) = tf(x,f), which is a weak solution of (1.1). In this 
N — *00 

section we show that this limit is in fact the unique entropy solution of (1.1) satisfying the 
entropy inequality (1.2). 

Theorem 4.1. Consider the spectral approximation ( 2.2 ) with parameters (e, m) which 
satisfy 

(4.1) e = e(JV) ~ Const. N ~ 20 , m = m(N) ~ Const. N 0 , 0 < /? < 1/4 . 

Assume that its solution, us(x,t) remains uniformly bounded, (3.17). Then un(z, t) con- 
verges boundedly a.e. to the unique entropy solution of the conservation law (1.1). 

Proof. Consider the righthand-side of (3.2) which consists of the sum of two terms, III + 
IV. We will show that this sum tends weakly to a negative measure and hence convergence 
to the entropy inequality (1.2) follows. 

As in (3.11) we write III = III\ 4- I IIi, where by Lemma 3.4 

(4.2) ||/J/i|| < Const. ||ttAr|| toi>(jt ^ • N 0 ~ l,i -> 0 , 

"toe 

and in view of our slightly strengthened parametrization (5.1) (compared with (3.18)), 

(4.3) ||/// 2 || t < Const.N 20 ~ l > 2 -* 0 , 0 < 0 < 1/4 . 

Consequently, the first term on the right of (3.2), III, tends weakly to zero, and we 
turn to deal with the second one which is given in (3.14) a s IV = IV \ + /V 2 . By Lemma 
3.5 we have 


(4.4) 


jTVx jj _ < Const. || us 




N~ 0 


0 , 


and hence the term IV \ also tends weakly to zero. Finally we are left with the term 
IV a = —e^£j*-{I - which we write as 


(4.5) 


d_ 

dx 1 


dun i 
dx 


d 2 u 


N 


dx 2 


= IV 11 + IV 12 + IV 23 



It follows from (3.20a), that —e{^££-) 2 tends weakly to a negative measure, 
(4.6) 


am [ JV„ = (5^.) 2 ] < 0 . 


w 

N—-oo 


Also, the pessimistic bound 

II d r d duN 


yields by Lemma 3.3, 


dx 




< £ • Ujv 




**(*.*) 


11 OX "t? (*,t) 


(4.7) ||/V M S < C« M ».||. w || t . (wi • N-* - 0 , 


and hence weak convergence to zero. We conclude with 

ao j%2 

(4.8) JV» = -eusF m = -gup m — s(un - = i ^231 + 1^232 • 


Here we have 


3 2 uw 

JV m = -£UP m -^f - 0 


(weakly), 


and by Lemma 3.3. 


3 2 tt> 


(4.10) 


||/V, 33 = - E (« w -u)P m -^f|| < 

loc 


< £• Utf - U 




m 


dx 2 l, A* <*, t) 


< Const. N • N 2 * 3 • 1 tzjy — u 




0 . 


In s ummar y, we have by (4.2)-(4.10) that the righthand-side tends in the sense 
distributions to a negative measure, while the (weak) limit on the left gives us 

Ifr’^r 3 ^ 0 ' 


as asserted. 
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5. Systems of Conservation Laws. 

It is instructive to compare between the spectral methods before and after spectral 
vanishing viscosity was added. Before viscosity was added we had, in (1.6), a coupled 
system of O.D.E’s in the Fourier space, which amounts to 

(5.1) j t u k (t) + l -ik £ 8p(()S,(l)=0 1*1 < N. 

p+q=k 

In this case the total quadratic entropy was conserved, (1.9), which is responsible to the 
divergence of the method. After viscosity was added in (2.2), the resulting system in the 
Fourier space reads 

(5-2a) + \ ik ^ “p(0“?(0 =0 |fc| < m . 

P+q=k 

(5-26) ^“*(0 + \ ik E “?(*)««(*) = -ek 2 u k (t) |fc| > m . 

P+q=k 

An increasing portion of the spectrum is treated here as in the diverging case (5.1). Yet, 
the added viscosity for the high Fourier modes in (5.2b) is responsible for the correct rate 
of entropy dissipation (3.6) which in turn implies convergence in the scalar case. In this 
section we show how to enforce similar entropy dissipation by spectral vanishing viscosity 
in systems of conservation laws. To this end we proceed as follows. 

Consider the conservative system 

(5.3) ^ tt (*,«) + ^[/(u(x,0)] =° > 

which is assummed to be equipped with an entropy function U(u), i.e., an convex function 
whose Hessian U uu symmetrizes the Jacobian matrix / u , e.g. [2], [11]. Using the entropy 
variables 

dU 

(5.4a) w = v(u) = -^-(u) , 

the conservative system (5.3) takes the equivalent symmetric form, consult [4], [8], [12] 


(5.46) 


§t u ( x >*) + = 0 , g(v) = /(«(„)) • 
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The Fourier approximation of (5.3) will be based on this formulation: together with addi- 
tional vanishing spectral viscosity we arrive at 


(5.5a) (*’*) + §^\ p N9{vN{x,t))] = 


{I - P m )-^v N {x,t) 


where vjy(x,t) is the projected vector of entropy variables 


(5.56) vs{x,t) = fV[v(ujv(s,t))] 

Multiply (5.5a) by v//(x,t) and integrate over the 2ir-period: taking into account orthog- 
onality we have 

d f w f* d 

— J U(utf(x,t))dx + J v N {x,t)—[g{v N {x y t)\dx 

/ 7f Q Q 

Vjv(x,t)^[(7 — P m ) ~q^ V n{x, t)]dx . 

The second integrand on the left is a perfect derivative of the associated entropy flux and 
hence its integral vanishes. Integration by parts on the right yields 

(5-6) U{u N {x,t))dx + e\\{I-P m )~v N {-^)\\\^ {t) =0 , 

which shows that entropy dissipates at the correct rate. In particular, arguing along the 
lines of Lemma 3.3 and using the strict convexity of U(u), we conclude that 

(5.7) \/£||^r^-|| a , , < Const. , em 2 ~ Const. 

This is analogous to the behavior of the viscous regularization for (5.3) - compare (3.20) 
in the scalar case. 


14 


0. Numerical Experiments. 

The Fourier method with spectral vanishing viscosity was applied to the periodic 
Burges’ equation (1.1) with u(x,t = 0) = stnx as initial data. The resulting O.D.E. 
system for the Fourier coefficients, see (5.2) 

(6.1a) ^Ufc(t) + it* ^ u p (t)u,(i) = -ek 2 Q(k)u k {t) , |*| < N , 

p+q=k 

<«*) «w-{! !*!;:: 

was integrated up to time t = 1.5, using the fourth-order Runge-Kutta method. 

The number of significant modes was chosen as a fraction of the total number, m — ON. 
The numerical experiments indicate, as expected, that the quality of the results is more 
sensitive to the dependence of s on m; further investigation is necessary in order to exhaust 
this point. In the following examples we have, em ~ 0.25, based on considerations of 
minimizing the total- variation of the numerical solution. With this choice of parameters, 
Figure 6.1 shows that the numerical solution converges strongly (but not uniformly) to 
the entropy solution of (1.1). This is in sharp constrast to the oscillatory behavior of the 
viscosity-free Fourier method in Figure 6.2, where Q(k) = 0*. Other parametarizations of 
e, quoted in Figures 6.3 and 6.4, demonstrate the sensitivity of the computed solution 
mentioned earlier. 

To improve the quality of these results, the proposed method (6.1a) was implemented 
with a spectral vanishing viscosity Q(k) which is smoothly varying between zero and one, 
say, for y < \k\ < m. Figures 6.5 and 6.6 show that this kind of viscosity prevents 
the propagation of the Gibbs phenomenan into the whole computational domain that was 
noticed earlier. This is analogous to the spectral recovery in shadowed regions between 
propagating linear discontinuties described in [7]. 


* In fact, a slight amount of dissipation was introduced in this case due to the time 


integrator. 
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The Fourier method with spectral vanishing viscosity ... and without spectral viscosity. 



Figure 6.1 


Figure 6.2 







The Fourier method with smooth spectral vanishing viscosity. 



Figure 6.5. 


Figure 6.6. 
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APPENDIX 

A. Failure of Convergence with Post-Processing. 

In practical applications the classical Fourier method is often coupled with certain 
smoothing procedures whose purpose is to gain spectral convergence that otherwise might 
be lost. In a typical case, the solution is post-processed via a convolution with a smoothing 
kernel Qn = Qw(s) 

N 

Qn * u N (x,t) = ^2 . 0 < < 9o = 1 > 9k = qk,N • 

k=-N 

In order to maintain spectral accuracy, the convolution with such smoothing kernel should 
be highly accurate with that of Dirac’s ^-distribution. We shall make a minimal assumption 
in this direction, requiring that for all square- integrable functions <p(x) we have 

(0.2 ) ||(«-Qw).*ip(*)|| l « ( , ) < ^2 °- 

k=-N ~~ C ’° 

Such smoothing procedure enables, for example, spectral recovery of solutions to linear 
hyperbolic problems in the presence of propagating initial discontinuities, e.g. [l], [7]. 

We will show that the smoothed version of the Fourier approximation to the nonlinear 
inviscid Burgers’ equation (1.1), 

d d 1 

(°.3) — u^(x,t) + ~q^\' 2 ^ >n ^ n * «*(x,t)) 2 ] = 0 > 

where the solution (rather than the flux) is convolved with kernels satisfying (a.2), does 
not converge to the entropy solution of (1.1). 

To this end we convolve (a.3) with Qn to find that u;^(i,t) = Qn * «v(x,t) satisfies 

d d 1 

(°-4) —w N (x,t) -I -Qn* 2 Pn^n{x^)] = 0 • 

Multiplying (a.4) by ujv(x,i) and integrating over the 27r-period we obtain 

/ UN ~W‘ dx ~~ J Qn * u n-^[^Pnu>n]<Ix , 
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or, in view of orthogonality 


“ J unQn * u N dx = - J w N ^-[^w%]dx = 0 . 


Thus with Qtf 2 * <p = ^2%=-n we ^ ave 

N 

(a.5) ||Qjr a *ttN(-.0||* = YL %\u k {t)\ 2 = Const. 

L (I) k--N 

This implies that our smoothed approximation t i>/v(x,t) = Qn * u N[x,t) converges to a 
weak limit wlim ws(x,t) = w(x, t). Now, suppose that w(x,t) is a weak solution of (1.1); 

N-+oo 

then this will lead us to the conservation of f g w 2 (x,t)dx which shows that w(x, t) is not 
the entropy solution (1.1). Indeed, if w(x,t) satisfies (1.1) then by (a.4), Qn * PNU>lf(x,t) 
tends weakly to W 2 (x,t) 


(o.6) w lim Qn * Pjvu;jy(x,t) = w 2 (x,t) . 

Pt—*oo 

It follows from (a.5) that ||Pjvu>jv( , ,t)|| i i( x ) is bounded, and together with (a.2) this 
implies that for all C^-test functions <p(x) we have 

(a. 7) j <p(x)(6-QN)*PNW 2 N{x,t)dx<Const.\\(6-QN)*<p{x)\\ L ' Xi -+0 . 

J X IS —>oo 

Adding this to (a.6), we conclude that /Vu/jy(x, £) and hence tt/jy(x,<) tend weakly to 
u7 2 (x,t). Consequently w(x,t ) is the strong limit of our smoothed approximation 


(a. 8) sliTOff-.ooQ n * u N {x,t) = w(x,t) . 

Finally, in view of (a.5) we can apply (a.2) to find that 

N 

(a.9) llQjv * u N (x,t) - Q]l 2 * ujv(x,t)|| tae(t) < 2 ^ |1 - ?*,jy| ■ qk,N • |«*(0I • 

From (a.8) and (a.9) it follows that W{x,t) is the strong limit of Q 1 ^ 2 * UN{x,t ), and the 
strong limi t of (a.5) tells us that f g w 2 (x,t)dx is conserved in time which completes our 
asserted contradiction. We summarize by stating 

Theorem A.l. The Fourier method (a.S) which employs any smoothing kernel satisfying 
(a.2), does not converge to the entropy solution of (l.l), (1-2). 
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